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ABSTRACT 

We derive a generalised van Cittert-Zernike (vC-Z) theorem for radio astronomy that is valid 
for partially polarized sources over an arbitrarily wide field-of-view (FoV). The classical vC- 
Z theorem is the theoretical foundation of radio astronomical interferometry, and its appli- 

I cation is the basis of interferometric imaging. Existing generalised vC-Z theorems in radio 

astronomy assume, however, either paraxiality (narrow FoV) or scalar (unpolarized) sources. 
Our theorem uses neither of these assumptions, which are seldom fulfilled in practice in ra- 

I dio astronomy, and treats the full electromagnetic field. To handle wide, partially polarized 

fields, we extend the two-dimensional electric field (Jones vector) formalism of the standard 
"Measurement Equation" of radio astronomical interferometry to the full three-dimensional 
formalism developed in optical coherence theory. The resulting vC-Z theorem enables all-sky 
imaging in a single telescope pointing, and imaging using not only standard dual-polarized in- 

I terferometers (that measure 2-D electric fields), but also electric tripoles and electromagnetic 

vector-sensor interferometers. We show that the standard 2-D Measurement Equation is easily 
obtained from our formalism in the case of dual-polarized antenna element interferometers. 

I We find, however, that such dual-polarized interferometers can have polarimetric aberrations 

at the edges of the FoV that are often correctable. Our theorem is particularly relevant to pro- 
posed and recently developed wide FoV interferometers such as LOFAR and SKA, for which 
direction-dependent effects will be important. 

Key words: telescopes; techniques: interferometric; techniques: polarimetric; instrumenta- 
tion: interferometers; instrumentation: polarimeters. 



Il995h which, in its original form, is a scalar theorjQ and only 
valid for narrow fields. Such restrictions are obviously unaccept- 
able in astronomy where polarization often provides crucial as- 
tronomical information, and sources are distributed on the celes- 
tial sphere, not necessarily limited to small patches. In radio as- 
tronomy, a polarimetric, or vector (i.e. non- scalar), extension o f 
the narrow-field vC-Z theorem was given bv lMorris et"ai] ( Il964l) . 
and more recently this wa s given a consistent mathematical foun- 
dation by iHamakej l lioool) in the form of the so-called Measure- 
ment Equation (M.E.) of radio astronomy. On the other hand, a 
wide-field extensi on of the scal ar vC-Z theorem for radio astron- 
omy was given bv lBrouwl ( fT97ll) . and more recently, imaging tech- 
niques for scalar, wide-fields have been an a ctive field of research 
JComwell & Perlev"l992!; ISault et al.|[T999l : IComwell et a"i]|2005l; 
[jVIcConnell et al. 2006). However, a generalised vC-Z theorem for 
radio astronomy that is both polarimetric and wide-field has not 
been derived. This may be because it has been assumed that such a 
theorem would be a trivial vector- or matrix-valued analogue of the 



^ By scalar theory we mean a theory which only considers unpolarized 
* E-mail: t.carozzi@physics.gIa.ac.uk radiation. 



1 INTRODUCTION 

Polarimetric wide-field imaging is a recent and important trend in 
radio astronomy. Many new and planned radio telescopes such as 
LOFAR, LWA, MWA, and SKA have a wide field-of-view (FoV) 
as a major feature. Technologically, this has been made possible 
by the development of phased dipole arrays and focal plane ar- 
rays, both of which have an inherently wider FoV than traditional 
single-pixel radio telescopes. The motivation for wide FoV polari- 
metric radio telescopes in astronomy is that they facilitate the study 
of polarized phenomena not restricted to narrow fields such as the 
large, highly-structured, pol arized features disc overed in recent po- 
larimetric galactic surveys jTavlor e t al. 2006). With LOFAR now 



producing its first all-sky images jLOFAR Teami.2007.) . a new era 
of polarimetric wide-field imaging is starting in radio astronomy. 

Despite this trend, a complete theory for wide-field, polari- 
metric astronomical interferometry is lacking. Ultimately, classical 
interferometry is based o n the far-zone form of the van Cittert- 
Zernike (vC-Z) theorem jThompson et alj|200ll ; iMandel & Wo3 
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wide-field, scalar theory jThompson et aljr200lh . As we will see, 
however, the final result is not that simple. 

The purpose of the present work is therefore to derive a vC-Z 
type relation by generalising the standard M.E. formalism to allow 
for arbitrarily wide fields. This is achieved by generalising the two- 
component Jones formalism to a three-component Wolf formalism 
l lWolj|l954h . In what follows, we will derive a vC-Z relation that 
is valid for the entire celestial sphere and is fully polarimetric. We 
will show that the standard M.E. can be recovered through a two- 
dimensional projection. We will also show that a dual-polarized 
interferometer of short electric dipoles (Hertzian dipoles) is inher- 
ently aberrated polarimetrically. We also extend our vC-Z relation 
to include not only the full electric field, but also the full magnetic 
field, and thereby establish the complete second-order statistical de- 
scription of the electromagnetic relation between source terms and 
interferometer response. 

Our wide-field vC-Z relations should be of relevance to 
the subject of direction-dependent effects in radio interferome tric 
imaging, which has attracted recent attention ( lBhatnagadl2008h . 



2 DERIVING A vC-Z RELATION IN 
THREE-COMPONENT FORMALISM 

The vC-Z theorem as used in astronomy is different from its orig- 
inal use in optical coherence. In astronomy one images sources 
on the celestial sphere based on localised measurements of their 
far-fields. This is possible because the vC-Z theorem provides an 
explicit relationship between the visibility measured directly by 
an interferometer and the brightness distribution of the sources 
jThompson et al.ir200ll . chap 14). This makes the vC-Z relation- 
ship the foundation of synthesis mapping and interferometric imag- 
ing. Despite its importance, the vC-Z relations in use in astronomy 
are all either based on the paraxial approximation (narrow FoV) or 
they take the source emissions to be scalar. These simplifications 
are questionable when the sources cover wide fields or are highly 
polarized, as is often the case in radio astronomy. 

Here we derive the full vector electro magnetic analogue to the 
wide-field, scalar vC-Z relation derived in [Thompson et all feOOll 
chap 14), and, as we will see, the final result is not simply a matrix- 
or vector-valued version of the scalar relation. We seek a relation- 
ship between the electromagnetic field coherence of a distribution 
of radio astronomical sources and the resulting electromagnetic 
field coherence at a radio astronomical interferometer. To simplify 
the discussion, we will use the term interferometer as a shorthand 
for radio astronomical interferometefl The treatment is intended 
for Earth-based interferometer observations, but it also has space- 
based radio astronomical interferometry in mind. 

Consider the problem illustrated in Fig.[T^). Radio emissions 
from far away sources, such as the point source S, are measured by 
an interferometer located in the domain M. We want to establish a 
relationship between the coherence of the electric field emanating 
from a source distribution and the coherence of the electric field in 
M. 

Let us first consider a single point source S located at R. The 
source is thus at a distance i? = |R| in the direction given by the 
unit vector (direction cosines) s — Tl/R with respect to the ori- 
gin O approximately at the centre of A4. The source is within the 

^ Most of the formalism also applies to single-pixel radio telescopes as a 
special case when the interferometer baseline length is zero. 



interferometers FoV which is centred on point C on the celes- 
tial sphere. The electric field from 5* is measured by the interfer- 
ometer at pairs of positions ri and r2 in A^. is assumed to 
be bounded, so the maximum distance between any two measure- 
ment points Dm = maxM jri — r2| (maximum baseline) is finite. 
Although the source emission may be broadband we split it into 
narrow, quasi-monochromatic, spectral bands and consider a typi- 
cal narrow (bandwidth much smaller than centre frequency) band 
centred on frequency u. The assumption that 5* is very far away, 
which quantitatively we take to mean that |R — r| 2> c/v for all 
r £ M, implies that the entire interferometer is in the far-field of 
the source. This means that the electric field E at point ri at time t 
is 

E(ri,s,t)=£fRi;s,i--^j — (1) 

where 

iii = |Ri| = In - i?s|, Ri=Ri/i?i (2) 

and £ £ C'"* is the complex electric field amplitude vector emitted 
by the source at s in the direction of ri . 

Furthermore, in the far-field, £ is approximately transverse so 

that 

Ri ■ f ^ 0. (3) 

If we further assume that the angular extent of M as seen from S 
is small, that is 

\ri\< Dm -tiR, (4) 

then the interferometer is in the Fraunhofer far-field of the sources. 
In this case we can use the approximations Ri ~ s and 1 /Ri « 
1/R to simplify the expression for the electric field in A4 to 

T^{r^,s,t)^£(s,t-ffj- (5) 

where 

s-£ ^0. (6) 

Here we have dropped the first argument of £ since we are as- 
suming that the angular variation of the sources emissions is 
small enough so that the Ri ~ s approximation implies that 
£ (Ri;s,t) ^ £ (s;s,t). 

A similar expression to equation Q for the field at r2 is ob- 
tained by replacing ri with r2 and Ri with i?2. The electric co- 
herence matrix (tensor) f can then be found by taking the outer 
product of the electric fields 

r,, ^{E,{ri,t)E*{r2,t)) (7) 

where E* denotes complex conjugation, ( ) denotes time aver- 
aging, and the subscripts i,j label Cartesian components x,y,z, 
which we will define more precisely later. We have written the elec- 
tric coherence matrix, equation ([7J, as a 3 x 3 complex matrix even 
though, for the single point source we are considering here, its rank 
is two and could therefore be expressed as a 2 x 2 complex matrix. 
We keep the electric coherence matrix as a 3 x 3 complex matrix 
since it is valid even when there are more than one point source. 

We now move to the case of a finite number of point sources. 
The total electric field measured at points ri and r2 is now the sum 
of the fields from sources in directions Sp for p — 1, 2, .... The 
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a) Position and visibility space b) Source direction, baseline and brightness space 



Figure 1. Coordinate systems and notations used in the text, a) A point source S located in look-direction s at a distance R with cross-section dA emits 
radiation, and its electric field E is measured at the pair of points ri and r2 located within the domain M. At the approximate centre of JV[ is the origin of 
the coordinate systems O. b) We use both a spherical and a Cartesian system to express the wide-field vC-Z. These systems are constructed with reference to 
C, the centre of the FoV If, and A an arbitrary direction orthogonal to C (which for Earth-based observations could be towards the North pole or zenith). The 
right-handed Cartesian system xyz is used for the source direction unit vector s = /x + my -f nz and the relative displacement vector D>, = nx + vy + wz,. 
The basis vector z points towards C and y points toward A. The angular part of the spherical coordinate system, {</>, 8}, has its zero point at C, and the pole 
j4 is at S = +7r/2. Relative to an arbitrary point on the celestial sphere given by s, the unit base vector 6 is directed along the arc from the point towards A, 
and ^ points in the anticlockwise direction as seen from A. These base vectors are used to specify the source elecfiic field e = e^(f> + eg6 and consequently 
also the two-dimensional brightness matrix. 



electric coherence matrix is therefore 

(n , t) (ra , t ) ) = ^ ^ (S. (r 1 , Sp , t) (ra , s„ t) ) 

p 1 

g-i2iri^s-(ri -r2)/c 

~ ('^* (^P' (sp. ^)) ^5 

P 

(8) 

for i,j = x,y, z. In going from the double to the single sum we 
used the usual vC-Z assumption that the sources are spatially inco- 
herent, that is, sources in different directions are statistically inde- 
pendent. 

Until now we have considered only discrete sources. We can 
make the transition to the more general continuum source distribu- 
tion by introducing the three-dimensional brightness matrix B'''' as 
a function of direction s in a continuous source distribution 

^ (f, (s„ t) E] (s„ t))^ \ B,<f (s) dA, (9) 

p J source 

where dA is the infinitesimal area of the source distribution. The 
superscript is to highlight the fact that B'"'^ is three-di mensional 
as op posed the usual two-dimensional brightness matrix ( lHamake3 
llOOOi) . The reason that the brightness matrix here is three- 
dimensional is simply because the full electric field amplitude £. is 
three-dimensional. However, due to equation l|6]l, not all the com- 



ponents of B^^ are arbitrary. In fact, we will show that it can be 
recast as one two-dimensional matrix. 

By making the replacement l|9j in equation l[8]l, we obtain 

(■ ... „-i2iri/s-(ri-r2)/c 

r(ri,r2,0)=/ B(^'(s)2 dA 

J source ^ 

= y Bf^)(s)e-''"''"("'i dJ7, (10) 

where AD. is an infinitesimal solid angle of the source distribution 
in T . In this expression we see that the dependence on the pairs of 
position is only relative, that is, the coherence matrix depends only 
on ri — r2 = D, and so we recast the expression in terms of the 
vector 

Da = -(ri - ra), 

c 

also known as the baseline vector measured in wavelengths (A = 

cjv). 

In practice, rather than use f directly, it is convenient and con- 
ventional to put the phase reference point at the centre of the FoV 
given by the direction sq . The result of the change in phase, 

V<^>(Da) = exp(i27rso ■ DA)r,(ri,r2), 

is the 3x3 generalisatio n of the standard 2x2 visibility matrix as 
defined by, for instance, iHamakeJ feOOOh . V''^^ is a general com- 
plex matrix except for Da = where it is Hermitian. It fulfill the 
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symmetry relation V*^'(-Da) = ^■^'^''(Da)^ where stands 
for Hermitian transpose. 

If we use V'^' in equation l llOt we arrive at 



,(3), 



B'=^>(s)e- 



i27r(s — so)-D;^ 



dn. 



(11) 



Equation i ll It is a matrix version of the wide-field, scalar vC-Z 
relation and so includes the partial polarization of the (non-scalar) 
source distribution. However, it is not very useful in this form since 
it does not automatically fulfill the constraint ((SJ. When applied to 
B'^\ this constraint becomes 







(12) 



for all directions s, where s is understood to be a column vector. 
The constraint can, however, easily be removed if we express £ 
in terms of spherical base vectors and set the radial component to 
zero. This leads us to introduce two coordinate systems: a Cartesian 
and a spherical. 

We use the angular, or tangential, basis set {0, 9} of a spher- 
ical polar coordinate system as the basis for the polarization of the 
transverse field of the source distribution on the celestial sphere, see 
Fig. lb). To simplify the results, the spherical coordinate system is 
taken relative to the phase reference position of the interferometer 
assumed to be at C. In other words, the zero point of the spheri- 
cal system, {</!>, 0} = {0, 0}, (intersection of the equator and cen- 
tral meridian), is taken to coincide with C. 6 is the angle from the 
equator (positive in the hemisphere with the pole A and negative 
in the other hemisphere), and <f) is the position angle from the cen- 
tral meridian around AO in the anticlockwise sense looking along 
AO. In lieu of any other reference directions, the orientation of the 
spherical system around C is arbitrary, but for Earth-based mea- 
surements A could be directed towards the North pole or zenith. 

One should note that {ci>, 0\ are consisten t with Ludwig's 
second definition as detailed in iLudwig' ('1973^ with the under- 
standing that the antenna boresight in Ludwig ( 1973) is here at 
C, that Ludwig's reference polarization unit vector i^^i is here 6, 
and Ludwig's cross polari zation unit vector icross is here </>. See 
jPiepmeier & Simorill2004l) for the use of Ludwig's third definition 
in a vC-Z relation. 

The Cartesian system, with base vectors {x, y, z}, is defined 
so that z is in the direction of C, and y is in the direction of the 
pole A. In terms of the Cartesian system, we can explicitly write 
the components of the vectors in equation i ll It as 



s = + my + nz = {I, m, n) 
so = z = (0,0,1)^, 
Da = wk + vy + wz = {u, v, w)"^ 



where 



(13) 
(14) 
(15) 



(16) 



and where the superscript ^ stands for vector transpose. All these 
vectors are unit vectors with real-valued components ranging be- 
tween — 1 and -1-1. These definitions of the uvw and Imn spaces 
are the same as the usual defini tions for Earth-based observations, 
see e.g., [Thompson eTZI j200lh . Also the matrices in equation i ll It 
are to be considered in what follows as being expressed in the 
Cartesian system. 

The relationships between the spherical and Cartesian systems 



base vectors are 



(nx — Iz) 



(17) 



(— Zmx + (1 — m^)y — mnz) . (18) 



Using these spherical and Cartesian systems we can express 
the three-dimensional transverse electric field as 



f = Te 



(19) 



where 



e{l, m) = 



e^{l, m) 
ee{l,m) 



is the Jones vector in spherical (rather than the usual Cartesian) 
components and 

.| / n —Im \ 

= 1 - for / 1 (20) 



—I —mn 



is the 3x2 transformation matrix between the components given 
by the equations l ll7t and l ll8t . Note that this transformation is pos- 
sible for all directions on the celestial sphere except for m — ±1, 
i.e., the poles of the spherical coordinate system. 

If one wishes to use to a polar spherical system in which the 
centre of the FoV is not on the equator, as it is here, but rather at 
some declination O, then one simply replaces T with T' defined in 
appendix|B] equation l lB6t . This assumes that the pole A is towards 
the Earth's North pole. 

It is easy to show that 

s'^Te = 0, 

so the transverse electric field expressed according to equation l |19t 
does indeed fulfill equation l|6ll. So if we use Te rather than £ in 
equation i ll It we would have an unconstrained vC-Z equation. We 
can introduce this replacement by rewriting B'^' using equation 
l fT9l l, so 

B(3) 



TBT" 



(21) 



B 



(22) 



where, suppressing the dependence on {I, m) 

(le^l^) {e^el) 
(ese*) (leel^) 

is the 2x2 brightness matrix, but in spherical rather than Cartesian 
coordinates. By this we mean that, for an arbitrary direction (I, m), 
B is locally equivalent to the usual paraxial brightness matrix in 
Cartesian coordinates. From its definition it easy to see that B is a 
Hermitian matrix. 

By using the Imn and uvw spaces, as spanned by the vectors 
s, and Da, we can write equation i ll It in a more explicit form. The 
exact form, though, depends on the extent of JT. If it is entirely in 
the hemisphere n > 0, then we write = ,7-+ and 

i27r[uZ + um+Tij(n — 1)] 



TBT 



T C 



■d/dm, (23) 



where we have used dil — dldm/\n\, and where the matrices de- 
pend implicitly on / and m. 

If, however, part of is in the n < hemisphere, then we 
must add to equation l |23t the contribution from this hemisphere 
given by the integral 

i27r[ui + T;m + ijj(n — 1)] 



// 



TB(/, m;n < 0)T 



T e 



■ d/dm, (24) 
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where T- is the subset of T in the n < hemisphere. 

The image horizon, n = 0, can also be included by 
reparametrising the integral in terms of (m, n) rather than {I, m) 
and using the replacement dQ — dmdn/j/|. The poles m — ±1 
can also be imaged, but one must then stipulate the orientation of 
the and 6 vectors at these singular points. 

Now, by extending to cover the entire n > hemisphere 
and extending JF+ to cover the entire n < hemisphere, the entire 
celestial sphere can be imaged in a single telescope pointing. An 
assumption here is that is a proper three-dimensional volume, 
or in other words, the baselines should be non-coplanar. If M is 
just a plane (coplanar baselines) then only one hemisphere can be 
mapped uniquely. 

Equation ([23} is our main result. It says that the full 3x3 
electric visibility matrix V'^' on the three-dimensional uvw space 
is given by the 2x2 brightness matrix B on the two-dimensional Im 
plane. The fact that this is a relationship between two matrices with 
different matrix dimensions is a fundamental feature of our vC-Z, 
and makes it clear that it is not just a matrix- valued generalisation 
of the wide-field, scalar vC-Z. Mathematically, this is ultimately 
due to the ranks of the fundamental matrices, of which we will 
speak more in section|4] 

The remaining vC-Z relationships that form a complete char- 
acterisation of the electromagnetic coherence response of a radio 
astronomical interferometer are given in appendix IaI 



3 WIDE-FIELD vC-Z RELATION AS A TRANSFORM 

In the previous section we derived a vC-Z relation, equation J23b . 
in which the visibility matrix is determined from the brightness ma- 
trix. It is well known that the original vC-Z theorem (far-zone form) 
for narrow-fields states that there is a two-dimensional Fourier 
transform relationship between visibility and brightness. In astro- 
nomical interferometry the transform aspect of this relationship 
is exploited to produce brightness images from measured visibil- 
ity. The wide-field vC-Z, equation l |23l l. is not a two-dimensional 
Fourier transform, but, as we will now show, it is still possible in- 
vert it and thereby establish a sort of generalised transform. 

First we should state that there are several ways of express- 
ing B in terms of V''^', even though the wide-field vC-Z relation, 
equation J23I ), is a one-to-one relationship in general. This is be- 
cause V'"^^ has redundancies, as one can expect considering the 
asymmetry in the respective matrix dimensions of B and V''^'. So 
although the rank of V'''' is three in general, it is overdetermined 
if B is known. We will first derive a solution that is valid for the 
entire celestial sphere based on the full V'^'. The case when only 
a projection of V*^"^^ is available will be discussed in section|5] 

Consider that we are given the full V'^' {u, v, w) and that we 
would like to solve equation ( 123b for B. B y extending the solution 
of the scalar problem described in lComwe ll & Perlev ( 1992) to the 
three-dimensional matrix relationship in equation i23\ we find that 
one approximate solution is 



B<^>(Z,m;rz > 0) = \/l - P - x 

M-M' 

(25) 

where A4 — A4' is the uvw space spanned by (r — r' ) / A for r G 
A4 and r' E Ai. From B'''', we can find a solution for the two- 



dimensional brightness matrix, 

B{l,m;n> 0) = T'^B'^'T. 



(26) 



An analogous expression applies for B'''-'(i,m;n < 0) but with 
the integration over n running from —1 to rather than to 1. 
For n = 0, an expression can be obtained by using the (m, n) 
parametrised vC-Z mentioned in the previous section. 
We can now state the generalised transform as 



V'''-'\u,v,w) B{l,m) 



(27) 



where ^ reads "is wide-field, polarimetric vC-Z related to". The 
relation from brightness matrix to visibility matrix is given (for 
n > 0) by equation {23}, and the relation from visibility matrix 
to brightness matrix is given by equations l l26b and ([25}. Note that 
the relation is not simply a two-dimensional Fourier transform 
as in the narrow-field case. Furthermore, the difference in the di- 
mensionality of V'^' and B make it clear that equation ([27} cannot 
simply be a matrix generalisation of the scalar vC-Z relation, as is 
sometimes assumed. 



4 USING STOKES PARAMETERS FOR WIDE FIELDS 

In practice it is common to use Stokes parameters to characterise 
the brightness and visibility matrices over narrow fields. Let us see 
how Stokes parameters can be applied to the wide-field vC-Z rela- 
tion, equation {23}. 

Let us first consider brightnesses. The difference between the 
two-dimensional brightness matrix B in equation ([23} and the usual 
two-dimensional brightness matrix in the paraxial approximation is 
that the former is defined on a spherical domain while the latter is 
defined on a Cartesian plane. Locally, for a source at some {I, m), 
the two brightness matrices can be made equal. Thus we can define 
the Stokes parameters in terms of the components of B with respect 
to a spherical basis {(jy,6} in analogy with Cartesian basis as 



/ 

Q 

U 

\ y J 



Bgg + B^^ \ 



(28) 



/ M ( 

Bee — B^^ 
Be<f, + B^e 
\ +i [Be^ — B^e) ) 

When the spherical system is aligned as it is in section |2] such that 
the intersection of its central meridian and its equator is located at 
the centre of the FoV and the pole A is towards the Earth's North 
pole, then in a sufficiently narrow field around the centre of the FoV 
the Stokes brightnesses (/, Q, ?7, are a pproximately equal to 
the Stokes parameters of the lAU faamak er & Bregman 1996). If 
one wishes to conform with lAU Stokes parameters over the entire 
celestial sphere, then one can use equation J28t with T replaced by 
T'(0), equation ( |B6b , where is the declination of the centre of 
the FoV. The correspondence is then that +0 is the lAU's +x, and 
+0 is the lAU's +y. 

Note however that the lAU basis set for polarimetry is not the 
same as the Cartesian basis set widely used in radio interferometry 
for defining source directions and baselines, which in this paper 
is denoted {x, y, z}. In adopting both these systems, therefore, a 
sacrifice must be made, and we have chosen to slightly modify the 
usual Pauli matrices that are used to relate the brightness matrix to 
the Stokes parameters. Explicitly, the brightness matrix in terms of 
the Stokes parameters in equation ([28} is in this paper 



B 



I-Q (7 + iV 
U -W I + Q 



(29) 
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As can be seen, the expansi on of this expression into unitary ma- 
trices, as in lHamake3 ( I2OO0I eq. 6), leads to matrices equivalent to 
Pauli matrices but with a change of overall sign for the matrices 
associated with Q and V. 

In terms of equation l |28l l, we can write the three-dimensional 
brightness matrix in equation l l23t as 



B 



(3) 



TBT' 



--IC1 + QCq + UCu + VC^ 



(30) 



where 



Co =- 



Cu =- 




It is easy to see that the d for i = I,Q,U,V depend only on 
{I, m, n), and that they become the three-dimensional analogs of 
the equivalent Pauli matrices in equation l l29t . 

Equation i30\ shows that for every look-direction {l,m), the 
three-dimensional brightness matrix has four degrees of freedom, 
here expressed as the four Stokes parameters. In other words, the 
four Stokes parameters completely characterise the partially polar- 
ized brightness also for wide fields, under the assumptions made in 
the derivation of equation l|23ll). 

Now let consider if Stokes visibilities can be extended to wide 
fields. A consequence of the three-dimensionality of V''^' in equa- 
tion l l23b is, however, that the Stokes parameters cannot in general 
fully characterise the electric visibility. This is because, in contrast 
to B'''', which has a rank of at most two due to the transversal- 
ity conditior[fl equation il2l . V'^' has no similar constraint. In- 
deed one can convince oneself of the full rank of V'^' by consid- 
ering two distinct point sources: the weighted sum of their B'"^^ 
matrices at some {u, v, w) point according to equation ([23}, will 
in general be rank three. As there are only four Stokes parameters, 
albeit complex-valued in the case of visibilities, they cannot fully 
parametrise a rank three matrix. 

Alternatively, rather than use the standard four Stokes parame- 
ters, one could also use com plexified versi ons of the nine, real, gen- 
eralised Stokes parameters jCarozzi et al. |[2000). These parameters 
are analogous to the standard Stokes parameters but can completely 
describe the coherence of the full three-dimensional electric field. 
In light of the discussion above, these generalised Stokes parame- 
ters are particularly suitable for parametrising V'"^' , but we will not 
discuss them here any further. 



^ That equation 112) implies that B^^^ has rank two comes from the rank- 
nullity theorem of linear algebra, since equation 112) implies that dimension 
of the null space is one and the matrix dimension of B'"'^ is three, so3— 1 = 
2 is the rank of B(^). 



5 GENERALISED MEASUREMENT EQUATIONS 

The vC-Z theorem is a basic, fundamental physical relationship is 
independent of technology. The measurement equation (M.E.) of 
radio astronomy, on the other hand, includes practical aspects of 
telescope measurements, in particular the instrumental response of 
the telescope. Usually it is a relationship between a 2 x 2 brightness 
matrix and a 2x2 cross-correlation matrix of the output- voltage 
of a dual-polarized interferometer. Since in the past such two- 
dimensional M.E. have been tacitly based on the paraxial approxi- 
mation valid only for narrow fields, it is important to verify that it 
can be recovered from the wide-field, polarimetric vC-Z, equation 
( 1231 ). for which the paraxial approximation is not used. Although it 
is possible to do this in a simple, straightforward way, we choose to 
do it in a more detailed way, introducing a formalism that extends 
the usual two-dimensional, electric field based model of radio as- 
tronomical antenna response. 



5.1 Electromagnetic antenna response model 

To obtain a 2-D M.E. we must first introduce a formalism for con- 
verting the full electric field to a voltage in the interferometer an- 
tenna. An electric field E at an antenna excites an open circuit volt- 
age V. Assuming linearity, these two quantities are related as 



V E 



(31) 



where L is the antenna effective length vector. In general it is a 
function of incidence direction, i.e. (/, m), but here we will only use 
ideal, Hertzian dipole antennas (short electric dipoles). These have 
the important property that their effective length does not depend 
on incidence direction, L is just a constant unit vector, and so it 
directly samples the component of the electric field along its length. 

If we have n co-located antennas that have no mutual cou- 
pling, their output voltages can be written in a matrix form 



V2 



( L 



\V„ J 



(1) 

X 
r(2) 



-(1) 
-(2) 



-(") 



'(1) 

'(2) 



(32) 



where the i-th row in the n x 3 matrix contains the components 
of antenna effective length vector L''\ The matrix of antenna ef- 
fective lengths, denoted L, has physical dim ension length and is a 
n-dimensional extension of the Q matrix in lHamaker et al.l j 19961 . 
eq. (2)), for which n = 2 and thus models dual-polarized an- 
tennas. When n = 3 and the antennas are linearly independent, 
then they sample the full three-dimensional electric field. Such 
an antenna system is called tri-polarized antenna in general, and 
tripole antenna if the three dipoles are approximately mutually or- 
thogonal. One can also include antennas that sample the magnetic 
field, and an arrangement of electric and magnetic antennas can 
be constructed so as to sample the full electromagnetic field at 
a point. Such antennas are called electroma gnetic vector-sensors 
dNehorai & Paldill 1 99 ll : [Bergman et a"l]|2005 !). In what follows, we 
will only be interested in dual- or tri-polarized antenna systems. 

In particular, let us consider a dual-polarized antenna that con- 
sists of two co-located, non-mutually coupled dipole antennas, one 
aligned along x and the other along y. The response of such a dual- 
polarized antenna is 



= L 



1 
1 



Ex 

Ey 

E, 



= L 



Ex 



(33) 
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where L is the antenna effective length. We have changed the sub- 
scripts on the voltages to reflect the right-hand side of the equation, 
in other words, for this dual-polarized antenna, each voltage com- 
ponent is directly proportional to a unique Cartesian component of 
the electric field regardless of the radiations incidence angle. Be- 
cause of this property, it can be regarded as an ideal dual-polarized 
polarimeter element. Associated with each dual-polarized antenna 
element is a two-dimensional plane in which the polarization is de- 
fined and measured. If this plane is the same for all of the elements 
in a dual-polarized interferometer or if all the planes are mutually 
parallel (common design goal for polarimetric interferometers) we 
will say that such an interferometer is plane-polarized, in analogy 
with plane-polarized waves. Note that a plane-polarized interfer- 
ometer is not necessarily a co-planar interferometer, and that a non- 
plane-polarized interferometer may be co-planar. If the plane of a 
plane-polarized interferometer is to be specified explicitly, we will 
say, e.g. in the case of equation l l33t . that it is a;j/-polarized. 

Although most existing polarimetric interferometers in radio 
astronomy are intended to be plane-polarized, it is possible to have 
more general antenna elements such as electromagnetic vector- 
sensor arrays or electric tripole arrays. An real-life example of the 
latter is the LOIS test station, see Bergman, Carozzi & Kar lsson. 
(2003); Bergman et al. ( 2005); Berg man & Carozzil ( l2008i) ; lThida 
(i2004,) ; ,Guthmann & Thide. C2005i) . The prime motivation for such 
a tri-polarized system is that it samples the full electric field in a 
single telescope pointing rather than just a projection. 

5.2 Recovering the 2-D M.E.: dual-polarized antennas and 
the paraxial hmit 

Now that we have introduced a model formalism for antenna re- 
sponse we can derive the 2-D M.E. for the special but important 
case of the xy-polarized (Hertzian dipole) interferometer, that is, 
the polarization plane is normal to the centre of the FoV. The out- 
put voltages from the xy-polarized elements located at points Pi 
and P2 are cross-correlated and the result can be expressed as the 
correlation matrix Rij(l, 2) =< Vi{l)Vj{2) > for i,j = x,y, 
where the arguments 1 and 2 refer to baseline points Pi and P2. L 
is matrix multiplied from the left and its transpose from the right 
with B'^-* , but since it does not depend on (/, m) in this case, it can 
be pulled out of the integral in ( 123 1 . and so the correlator output can 
be written 



R 



{^y)^}{'■)) 



where 



1 
1 



(34) 



is the a;j/-polarized antenna elements, effective length matrix. As 
one can see, the effect of L*-^^-* is equivalent to projecting the field 
vectors into the xy-plane and multiplying by L. In terms of the 
brightness matrix the correlator output is 

R = ILpx 



rp —i27T\ul-\-vm-\-w(\/l — l'^—in'^ — l]\ 

f j t(..)b {j^-y)^ ^ — . dWm 



(35) 



where we have introduced the xy-projected transformation matrix 



-p(a:t/) _ _;_^_{xy)-J■ _ ^ 



n —Im 



to simplify the final result, which is now clearly two-dimensional. 
Equation ([35} is a wide-field M.E. but with the novel Jones ma- 
trix T'^"' that physically represents a projection of the three- 
dimensional electric field vector onto the xy-plane. 

In the narrow FoV limit \/P + m? <^ n, so we can approx- 
imate T^^^' in equation i35i up to first order in I and m by the 
two-dimensional unity matrix, and so 



R = \L\ 



il 



Be"' 



27r(uZ + um) 



dldm for ^/P + m? <^ 



(36) 

This is the basic M .E. of astronomical interferometry, see 
[Thompson et al.' ( 1200 iL eq. (14.7)). Thus we have shown that the 
wide-field, polarimetric vC-Z, equation i23i indeed reduces to the 
usual two-dimensional, Jones vector based M.E. in the paraxial 
limit. The result, equation i36h depended on the particular the 
spherical coordinate system used as default in this paper. Only this 
particular choice reduces directly to the standard 2-D M.E. in the 
paraxial limit. 

To see what the dual-polarized M.E. equation ([36} misses by 
not including the third dimension along z, let us go back to equation 
( 1231 ) and let assume the paraxial approximation, y/P + m? <C n. 
In this case 



(37) 



where we have kept only terms of first order in I and m. V^^ in 
equation ([23} is identical to Rij/\L\^ in equation ([36} for i, j — 
X, y, but for the z components. 




V. 



and 



(3) 



B 



VS' ^ 0. 



i = x,y 
(38) 

(39) 



So the zz component of the three-dimensional visibility matrix 
does not provide anything, but the xz (and zx) and the yz (and zy) 
components do. Thus even for a narrow FoV, a dual-polarized inter- 
ferometer does not measure the full set of generally non-zero elec- 
tric visibilities. Note that if the antenna array is not exactly plane- 
polarized, these additional visibility components will contribute to 
the output-voltage of such an array. 

This leads to the important question of how serious the loss 
of visibility information is in a dual-polarized interferometer. More 
specifically, we ask whether a plane-polarized interferometer, given 
by some 2x3 matrix L, can in general perform full polarimetry 
of an arbitrary source distribution. The answer is that the source 
brightness matrix in some direction can be determined fully only if 

det (LT) / 

since in this case LT is invertible. For the special but important case 
L = L'^^' this condition is equivalent to 

Thus, an sy-polarized interferometer can recover the full polarime- 
try except on the great circle orthogonal to C, which we may call 
the imaging horizon. So, in wide-field imaging with a;j/-polarized 
interferometers the image horizon cannot be measured with a single 
telescope pointing. Under noise-free conditions, this would pose 
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little problem, but if we add in the effects of noise then the great 
circle of directions for which full polarimetry is not feasible broad- 
ens as a function of the signal-to-noise ratio. 

Although the discussion above was mainly focused on plane- 
polarized interferometers, the three-dimensional formalism devel- 
oped here can also be applied to more general dual-polarized in- 
terferometers. In particular it can model the situation when the po- 
larization planes of dual-polarized antennas in an array are not all 
parallel. Since several large arrays of dual-polarized antenna based 
interferometers are currently being planned for, this would address 
the important question of whether such arrays should strive to be 
plane-polarized or whether they should purposely not be plane- 
polarized to minimize the inversion problems mentioned above. 



5.3 Polarization aberration in xy-polarized, Hertzian dipole 
interferometers 

We now show that the M.E. for the a;?/-polarized, Hertzian dipole 
interferometer, equation l |35t , exhibits distortions that depend on 
the look-direction, that is, the i mages based on these bright nesses 
contain polarization aberrations jMcGuire & ChiDmanlll990l) . This 
agrees with the general understanding in observational radio as- 
tronomy that the polarimetry of a telescope is worse off-axis than 
on-axis. 

Say we have measured the 2x2 correlation matrix R{u, v, w). 
We cannot use the formal solution ( 125 1 directly since it is for the 
3x3 visibility, and it is not clear how to obtain the remaining, un- 
measured z components from the x and y components of R. On 
the other hand, for the scalar case, the formal solution for produc- 
ing a synthesized image is the following scalar, wide-field imaging 
equation: 



(40) 

where / — 1(1, m) is the scalar brightn ess and V = V(u, y, w) is 
the scalar visibility, see equation (13) in IComwell & Perlevl l ll992h 
who call it the 3D method of synthesis imaging. Let us apply ( 140) 
to each component of R as if it were a scalar, thus creating a matrix 
analogue of the scalar, wide-field imaging equation. The resulting 
brightness matrix B'^^' (polarized image) is 



\L\' 



I I I I w~» 27vi\ul4-vm-\-'w(n — l)] i ^ 

■JJJjRe )ld«d 



vdwdn. 
(41) 



However, this is not the true brightness matrix because we see from 
equation ( I35t that 



B 



It follows that B'^''' is actually the projection of the three- 
dimensional brightness matrix into the xy-plane. The Stokes xy- 
projected brightnesses Sp are based directly on B'^"^ in analogy 
with equation ( I28t . that is. 



Sp{l,m) 



I. \ 
Qp 
Up 
Vp J 



/ 




1 r>(^v) 






^yy 




J-^xy 


~r J->yx 


\ 


i {b17^ 


r>(xy) 
J^yx 



(43) 



The relationship between these Stokes brightnesses and the true 
Stokes brightnesses S, which are based on B can found by recasting 



equation i42l as 



Sp = MS, 



(44) 



where 



/ i(l + n^) 



M 



(1 

Hi' 



2{l-m^) 
_ (l + m^)i^ 



V 



-Im 




2(1- 

-Im 




m2) 



n 




/ 
(45) 

M is a Mueller matrix that quantifies the distortion of true Stokes 
brightnesses based on the imaging equation ( I41l >. 

For a very narrow FoV, {I, m) m (0, 0) and M is approxi- 
mately unity. In general, however, M is not the unit matrix with 
the effect that the perceived Stokes vector is a distortion of the true 
Stokes vector. Thus, without further processing a plane-polarized 
interferometer of short dipoles will exhibit polarization aberrations 
over wide fields. By contrast, a tripole array interferometer is, at 
least in theory, polarimetrically aberration-free over a wide field, 
since the scalar wide-field imaging equation J40I ) applied to the 
components of its visibility matrix as if they were scalars gives 
B'^\ which can be interpreted as the exact B in Cartesian (rather 
than spherical) components. 

Examples of these wide-field distortions are displayed in Fig. 
|2] It shows aberration effects for source distributions that are con- 
stant over the entire hemisphere, by which we mean that B(7, m) = 
B, so the brightnesses do not explicitly vary with directiorQ Fortu- 
nately, these distortions can be compensated for in the image plane 
since M is invertible and well-conditioned as long as + is not 
close to one. 

Another aspect of equation ( I44t that is important, is that it 
shows, not surprisingly, that partially polarized, wide fields cannot 
be treated as scalar, wide fields. In fact, even for the supposedly 
'scalar' case, that is, when we only consider the scalar visibility 

V{u, V, w) 
Q = U = V 



2|jp-Tr(R) of an unpolarized source distribution, so 
0, equations ( 1441 and ( BTt imply that 



V{u,v,'w) ^ 1(1, m)- 



2-r 



(46) 



where ^ stands for two-dimensional Fourier transform. This 
should be compared with ^^(it, w, ui) I {I, m)/\/l — P — m? 
for th e isotropic, scalar antenna case (see, e.g., equation (1) in 
IComw ell & Perlev ( 1992)). This expression differs even for narrow 
fields since, to lowest, non-vanishing order in I and m, equation 
( |46l l becomes approximately 



V{u,v,w) ^ 1(1, m) 



1-1(1^+ my 



(47) 



while V ^ J(l — 1/2(Z^ + m^)) for the isotropic, scalar antenna 
case. Thus short dipoles are less aberrated than isotropic, scalar 
antennas. 

The conclusion here, that scalar theory is not sufficient for 
the description of the general vC-Z relations, agr ees with thos e 
detailed in the field of op tical coherence, see e.g. ICartetl ( Il980l) : 
ISaastamoinen et al] (|2003'). Also, equation i46l corresponds to an 
analogous equation in Saastamoinen et al. (2003). Consider an un- 
polarized point source at {l,m,n) — {I' ,m' ,n'), so 1(1, m) = 



This does not mean that the distributions are necessarily isotropic since 
they can still vary through their reference to the (direction-dependent) basis 



{0, 0}. 
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Figure 2. Distortion of various polaiized source distributions across tlie 
hemispliere for an xy-polarized interferometer All distributions are such 
that the Stokes brightnesses are constant, that is, they do not vary explic- 
itly with direction although they may vary implicitly due to variation of the 
reference system for the Stokes brightnesses, and cf>, with direction. The 
plots are of the polarization ellipses corresponding to the normalized Stokes 
parameters {Qp/Ip, Up/Ip, Vp/Ip). where {Ip, Qp, Up, Vp) are the xy- 
projected Stokes brightnesses as a function of the direction [l, m). The plots 
clearly show polarization aberration, that is, a direction dependent distor- 
tion in the observed polarization. The source distributions state of polariza- 
tion (shown in red) can be seen at the centre of the FoV, (I, m) = (0, 0), 
where there is no distortion. The upper-left panel shows S = (1, 0, 0, 0), 
i.e. completely unpolarized radiation, the upper-right panel shows circularly 
polarized radiation, the lower-left shows radiation linearly polarized along 
0, S = (1,1, 0, 0), and the lower-right panel is for radiation linearly polar- 
ized along 0, S = (1, -1, 0, 0) . 



IS{1 — l',m ~ m')j7i|, and express n' as cosa (where a is the 
angle between the point source po sition and C, and is equivalent 
to 6 in lSaastamoinen et alj ( |2003|) ), then the intensity of the un- 
polarized point source, as measured by single-pixel telescope, is 
V (0, 0, 0) = ^Jfl+ cos^g ). This corresponds to equation (47) 
in Saastamoinen et alj ( l2003h . 

Although the results will be different for other types of an- 
tennas, the Hertzian dipole is an important special case as it is the 
simplest polarimetric antenna and they are directly proportional to 
the Cartesian coordinates of the electric visibility matrix. 

The ultimate reason for the aberration is of course that the 
incident, transverse field is being projected onto the polarization 
plane of the Hertzian dipoles thereby distorting the field. Indeed 
this projection is identical to the orthographic projection of a hemi- 
sphere in cartography. In comparison to other effects encountered 
in instrumental calibration, such as beam-shape, these effects may 
seem small. They are, however, important in that they determine the 
ultimate limits of polarimetry since these effects are of a intrinsic, 
geometric nature. 

Inspection of M reveals that around {I, m) = (0, 0) the aber- 
rations are O (Z^ , ) . In the case of the aperture a rray of the SKA, 
the preliminary specification dSchilizzi et al.l2007h has a mean FoV 



of 250 square degrees, so the aberration error is in the order of -2%, 
or - 16 dB, at the edge of the FoV. This should be compared with the 
requirements from the key scientific projects that in some cases call 
for -30 dB in polarization purity over wide fields. Thus, these polar- 
ization aberrations will need to be considered in wide-field imaging 
with SKA. 

It is not difficult to show that these polarization aberrations 
occurring at the edges of wide-field images also occur for off-axis 
imaging with fixed mount telescopes. This occurs, for instance, in 
phased arrays of crossed dipoles such as LOFAR and the low fre- 
quency part of SKA, where imaging at scan angles away from bore- 
sight (zenith) will be achieved by electronically steering the beam- 
form. Since the crossed dipoles are fixed to the ground and are not 
on a mechanically rotating mount, a situation geometrically analo- 
gous to the wide-field imaging consider above occurs. This leads to 
similar aberrations in polarization for short electric dipole arrays. 
Scan angles of 45° have been proposed for which the aberrations 
will of the order of -3 dB. 



6 CONCLUSION 

We have derived the full set of electromagnetic vC-Z relations, 
which are the basis of radio astronomical interferometry, without 
invoking the paraxial approximation. These relations allow all-sky 
imaging in a single telescope pointing. We have achieved these re- 
lations by generalising the usual 2-D Cartesian Jones vector based 
M.E. to a 3-D Wolf coherence matrix formulation on the (celes- 
tial) sphere. The derived wide-field vC-Z relations are not simply 
trivial matrix (or vector) analogues of the wide-field, scalar vC-Z, 
they exhibit direction dependent polarimetric effects. Indeed even 
in the scalar limit (that is, unpolarized radiation), our M.E. is not the 
same as the M.E. derived from scalar theory, in the case of Hertzian 
dipoles. Furthermore, we found that, for an arbitrary wide-field of 
sources, the electric visibilities generally have nine complex com- 
ponents for an arbitrary baseline. This implies that the standard 
Stokes (Cartesian) visibilities do not provide a full description of 
electric coherence for wide fields. We have also shown that our 
vC-Z relation (for the electric field) reduces to the standard 2-D 
M.E. after a 2-D projection. We showed that a consequence of this 
projection is that plane-polarized, Hertzian dipole interferometers 
are aberrated polarimetrically. Fortunately, these aberrations can be 
corrected for in the image plane for sources sufficiently far away 
from the plane of the dual-polarized antenna elements. 

Besides its use in the derivation of the wide-field vC-Z rela- 
tion, we believe that the 3-D Wolf formalism can be useful in con- 
structing more general 3-D M.E. in cases requiring the full set of 
electric components. As examples, we mention the modelling of tri- 
polarized element arrays, the modelling of dual-polarized element 
arrays that are not strictly plane-polarized (due to manufacturing 
errors or the Earth's curvature), and the modelling of propagation 
that is not along the line-of-sight (due to refraction or diffraction in 
the ionosphere, e.g.). 
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APPENDIX A: 
RELATIONS 



FULL ELECTROMAGNETIC vC-Z 



In the previous sections we considered only the electric field and 
its auto-correlation, but now we look at the full electromagnetic 
field. One motivation to do this is that, especially for low radio 
frequencies, it is possible to sample both the electric and the mag- 
neti c fields and thereby measure e l ectrom a gnetic coherence fully , 
see 



: netas ana thereby measure e l ectrom a gnetic coherence tulty , 
Bergman. Carozzi & KarlssonI ( l2003h . iBergman et alj ( l2005h 



iBergman & Carozzil ( l2008l> . Such sensors have been deployed 
in some radio interferometers and will possibly provide unique 
and novel astronomical measurements. For completeness then, we 
present the rest of the electromagnetic correlations analogous to the 
electric vC-Z relation in equation l l27t . 

We found previously that the 3x3 electric visibility matrix 
was related to the 2x2 electric brightness matrix as 

V^'"'\u,v,w) ^ B^"'\l,m) (Al) 

where denotes the vC-Z relationship given by equation i23l . 
Note that we have now changed the notation of V'^' to v'^^' 
and B to B'^'^'^' to indicate that these quantities represent auto- 
correlations of visibility electric field and the brightness Jones 
vectors, respectively. A full electromagnetic vC-Z relationship 
between brightnesses and visibilities requires also the auto- 
correlations of the magnetic field and the cross-correlation between 
the electric and magnetic fields, see Mandel & W olf (1995, chap. 
6). We define the electromagnetic visibilities for baseline Da with 
respect to the phase reference direction so as 

Vif^^(D,) = (S,(r)S*(r- 



V, 



V 



(HE), 



D)) exp(i27rso • 
' D)) exp(i27rso 



V, 



Da) 
Da) 
Da) 
Da) 



(A2) 
(A3) 
(A4) 
(A5) 



D)) exp(i27rso 

. '(Da) = {H,{r)H;{r - D)) exp(i27rso 

for i,j — x,y, z, where H(r) is the magnetic field at r. The elec 
tromagnetic brightnesses in direction s are defined as 

B^f (s) = (e,(s)e*(s)) 
Sif'(s) = <ft.(s)e*(s)) 



(A6) 
(A7) 
(A8) 
(A9) 



for i,j = (j),0, where h(s) is the magnetic field from the source 
distribution in direction s. 

The magnetic field associated with the visibilities can be de- 
rived by applying Faraday's law, 

H(r,t) = -i^VxE(M), 

where Zo is the impedance of free space, to the electric field used 
in derivation in section[2] The result is that the magnetic analogues 
of the electric vC-Z expressions involve a magnetic Jones vector 
h — {hrf,, hg)'^ that is related to the electric Jones vector through 

1 



-s X (Te) = TFe = Th 



where 



or in other words 



1 

-1 



(AlO) 



Wide-field, polarimetric van Cittert-Zernike theorem 1 1 



This says that the magnetic Jones vector is directly determined from 
the electric Jones vector, and so no other independent electromag- 
netic source coherence statistics exist (in the far-field zone) other 
than the 2x2 electric brightness matrix. 

Repeating the derivation in section[2]but for the magnetic field 
we find we can write the rest of the vC-Z relations as 



g(eh) ^ g(ee)pT 



,{HH) 



(he) 



(hh) 



FB 
FB 



(ee) 



(ee)pT 



(All) 
(A12) 
(A13) 



Thus one can see F as a Jones-like matrix that switches between 
electric and magnetic coherencies. 

These relations provide a complete description of the second- 
order coherence of the electromagnetic radiation. So, for instance, 
one can easily compute the Poynting visibility vector using the 
above relations, 



[E (r) X H* (r - D)]) 



(A14) 



where / = /(/, m) — |eep + le^l^. This says that the power flux 
visibility vector is vC-Z related to the Stokes brightness propagat- 
ing from sources. 



APPENDIX B: COORDINATE TRANSFORMATIONS 

In this paper we have used the standard definitions of the uvw- 
and imn-spaces, see Thompson et al. ( 2001), and we have used a 
spherical basis {6, 0}. Both these coordinate systems are relative 
to the centre of the FoV. Often one wants to transform to some 
other system, and som e details can be found in the standard texts 
('Thompson et al."200l'). However, what is not usually done explic- 
itly is the transformation of the non-scalar brightnesses and vis- 
ibilities. We will now show how to rotate the matrix brightnesses 
and visibilities found in the polarimetric, wide-field vC-Z, equation 
([23}. 

Equation J23I ) uses the Cartesian {x, y, z} system as a basis 
forV'^\s,DA and the rows of T. It uses the spherical base vectors 
{0, 4>} as a basis for B and the columns of T. Under a rotation 
given by the 3x3 orthogonal matrix Q, the vectors and matrices in 
the Cartesian system transform in the usual manner: 



s' = Qs, 
D'a = QDa, 
V'(^^(s') =QV(^)(s')Q' 



(Bl) 
(B2) 

(B3) 



Note that a transformation analogous to equation ( IB3t for the 2x2 
matrix R, discussed in section [5!2l does not exist since R lacks one 
of the dimensions necessary for a general coordinate transforma- 
tion. Only the correlator output from a tri-polarized antenna array 
can be fully transformed in general. 
The brightness matrix rotates as 

r(s')QT(s)B(s)T^(s)QTT(s'), 



B'(s') 



(B4) 



usmg 



where s on the left-hand side should be converted into s 
equation dBlt . 

Actually, Q need not be the same in equations JB3I > and JB4I (. 
since the spherical and Cartesian systems can be rotated separately. 
This can be used to change the relative alignment (rotation) of the 
spherical system relative the Cartesian system. The net result is a 



change in matrix T that relates the spherical components to the 
Cartesian components. Consider the special case when the spheri- 
cal system is rotated in the positive sense around the x-axis through 
angle G relative the Cartesian system. The effect on the vC-Z rela- 
tions is that T is replaced by T', where 



(B5) 



— (m cos — 71 sin O)-^ 

(n cos Q + m sin O 
— ZsinG (1 — m'^) cos + mnsin O 

—I cos Q ~mn cos + (n^ — 1) sin O 



-Im cos O + in sin O 

,2\ 



(B6) 

Obviously, for O = we obtain T' = T, which is the matrix used 
in most of this paper. 

At the field centre, / = m = 0, so 



(B7) 



for all O 7^ 7r/2. However, the x and y components of its first 
derivatives are zero only for = 0, that is 




9T', 



dl 



= 0, 



aTL 



dm 



, t — x,y 
= 0, for <^ . ' ^ 



only for = 0. Thus, the special case we have used in this paper, 
T' = T, can be said to possess a projection (xy) that is locally 
flat at the field centre. This is a reason for choosing the spherical 
system with O = as a default, since all other cases would lead to 
a 2-D M.E. l l36t with additional first-order terms that account for 
the coordinate system curvature within the FoV. 

The matrix T' can be used to adapt the vC-Z relations given 
in this paper, such as equation l l23t , to standard celestial coordinate 
systems such as the equatorial system or the azimuth-elevation sys- 
tem. In essence, assuming that the pole A is towards the Earth's 
North pole (in the case of equatorial coordinates) or zenith (in the 
case of az-el coordinates), one simply performs T — * T' and then 
interprets O as either the declination (equatorial case) or the ele- 
vation (az-el case) of the centre of the FoV. Note however that the 
Cartesian coordinates, in which the uvw and Imn spaces are ex- 
pressed, would still need to transformed for complete agreement 
with these standard celestial s ystems, but for this com mon task we 
refer to standard texts such as [Thompson et al] ( 1200 ih . 



